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ABSTRACT 

The purpose of this study is to investigate compressibility effects on the turbulence in homo- 
geneous shear flow. We find that the growth of the turbulent kinetic energy decreases with 
increasing Mach number - a phenomenon which is similar to the reduction of turbulent ve- 
locity intensities observed in experiments on supersonic free shear layers. An examination of 
the turbulent energy budget shows that both the compressible dissipation and the pressure- 
dilatation contribute to the decrease in the growth of kinetic energy. The pressure-dilatation 
is predominantly negative in homogeneous shear flow, in contrast to its predominantly pos- 
itive behavior in isotropic turbulence. The different signs of the pressure-dilatation are 
explained by theoretical consideration of the equations for the pressure variance and density 
variance. We obtained previously the following results for isotropic turbulence; first, the 
normalized compressible dissipation is of 0(A/ t 2 ), and second, there is approximate equipar- 
tition between the kinetic and potential energies associated with the fluctuating compressible 
mode. Both these results have now been substantiated in the case of homogeneous shear. 
The dilatation field is significantly more skewed and intermittent than the vorticity field. 
Strong compressions seem to be more likely than strong expansions. 


1 This research was supported by the National Aeronautics and Space Administration under NASA Con- 

tract No. NAS 1-18605 while the authors were in residence at the Institute for Computer Applications in 
Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665. 
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1 Introduction 


Homogeneous shear flow refers to the problem of spatially homogeneous turbulence sustained 
by a parallel mean velocity field u = (£ 2 : 2 , 0,0) with a constant shear rate S. Such a flow 
is perhaps the simplest idealization of turbulent shear flow where there are no boundary 
effects, and where the given mean flow is unaffected by the Reynolds stresses. Nevertheless, 
the crucial mechanisms of sustenance of turbulent fluctuations by a mean velocity gradient, 
and the energy cascade down to the small scales of motion are both present in this flow. 

The homogeneous shear flow problem has been studied experimentally by Champagne, 
Harris and Corrsin (1970), Harris, Graham and Corrsin (1977), and Tavoularis and Corrsin 
(1981) among others. In these low-speed experiments the statistical properties of the flow 
do not vary spatially in the transverse (x 2 ,X 3 ) plane but they evolve in the streamwise x\ 
direction. In theory, the streamwise inhomogeneity can be removed by writing the equations 
of motion in a reference frame moving with the mean flow IT. In the moving reference frame, 
the one-point moments satisfy pure evolution equations in time clearly illustrating that ho- 
mogeneous turbulence is fundamentally an initial value problem. Rogallo (1981) and Rogers 
and Moin (1987) have investigated the incompressible homogeneous shear problem at great 
depth through direct numerical simulations. These simulations, albeit at low turbulence 
Reynolds numbers, have provided turbulence statistics which are in good agreement with 
experiments performed at relatively higher Reynolds numbers. Furthermore, since the sim- 
ulations provide global instantaneous fields, the turbulence can be studied in much greater 
detail than in physical experiments. 

Recently there has been a spurt of activity in the direct numerical simulation (DNS) of 
three-dimensional compressible turbulence. Decaying isotropic turbulence has been studied 
by Passot (1987), Erlebacher et al. (1990), Sarkar et al. (1989), and Lee, Lele and Moin 
(1990). The simulations of Erlebacher et al. (1990) identified different transient regimes 
including a regime with weak shocks, and also showed that a velocity field which is initially 
solenoidal can develop a significant dilatational component at later times. Sarkar et al. 
(1989) investigated the statistical moments associated with the compressible mode in their 
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simulations, and determined a quasi-equilibrium in these statistics for moderate turbulent 
Mach numbers which was then used to model various dilatational correlations. Lee, Lele 
and Moin (1990) studied eddy shocklets which developed in their simulations when the 
initial turbulent Mach number was sufficiently high ( M t > 0.6). Kida and Orszag (1990) 
primarily studied power spectra, and energy transfer mechanisms between the solenoidal and 
dilatational components of the velocity in their simulations of forced isotropic turbulence. 

Physical experiments have not been and perhaps cannot be performed for homogeneous 
shear flows at flow speeds which are sufficiently high to introduce compressible effects on the 
turbulence. However, direct numerical simulation of this problem could provide meaningful 
data, especially since DNS for the incompressible problem has been successful in giving 
realistic flow fields. The compressible problem was considered by Feiereisen et al. (1982) 
who performed relatively low resolution 64 3 simulations and concluded that compressibility 
effects are small. Recently Blaisdell (1990) has also considered compressible shear flow. 
We have performed both 96 3 and 128 3 simulations which have allowed us to obtain some 
interesting new results regarding the influence of compressibility on the turbulence. In 
contrast to the results of Feiereisen et al., our simulations which start with incompressible 
initial data develop significant rms levels of dilatational velocity and density. We find that 
the growth rate of the kinetic energy decreases with increasing Mach number as well as 
increasing rms density fluctuations and show that the compressible dissipation and pressure- 
dilatation contribute to this effect. Apart from rms levels of the fluctuating variables, we 
have also examined their probability density functions (pdf) and higher order moments. The 
dilatational field has different skewness and flatness characteristics relative to the vorticity 
field. In particular, the dilatational field has a significant negative skewness and is more 
intermittent than the vorticity field. 

2 Simulation Method 

The compressible Navier-Stokes equations are written in a frame of reference moving with the 
mean flow Uj. This transformation, which was introduced by Rogallo (1981) for incompress- 
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ible homogeneous shear, removes the explicit dependence on u(x2) in the exact equations for 
the fluctuating velocity, thus allowing the imposition of periodic boundary conditions in the 
x 2 direction. The relation between x* and the lab frame x* is 

Xj = X\ — Six 2 , = X2 , x 3 = X3 

Here S denotes the constant shear rate In the transformed frame x^, the compressible 
Navier- Stokes equations take the following form 

dtp + (pui),i - St(pu 2 '), 1 = 0 (1) 

d t {pu x ') + ( pu/ui )j = — P,i + T ij,j — Spu 2 Sil 

+ St{pu 2 'u x , ') t i + Stp t i 6 i 2 Str ^ 2,1 (2) 

d t p + u/pj + 7 pu'jj = Stu 2 'p tl + jStpu' 2l 

+ (7 - 1M Tjj - 2 StT l2 + S 2 t 2 T,n] + ^ (3) 

p = pRT (4) 

where $ = TijUij is the dissipation function, Ui the fluctuating velocity, p the instanta- 
neous density, p the pressure, T the temperature, R the gas constant, and k the thermal 
conductivity. The viscous stress is 

2 

Tij = + u jti ) - -p(u ktk )Sij 

where p, is the molecular viscosity which is taken to be constant. All the derivatives in the 
above system are evaluated with respect to the transformed coordinates x t . 

Since Eqs. (l)-(4) do not have any explicit dependence on the spatial coordinates x*, and 
because the homogeneous shear flow problem, by definition, does not have any boundary 
effects, periodic boundary conditions are allowable on all the faces of the computational box. 
Of course, in order to obtain realistic turbulence fields it is necessary that the length of the 
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computational domain be much larger than the integral length scale of the turbulence. Spec- 
tral accuracy is obtained by using a Fourier collocation method for the spatial discretization 
of the governing equations. FFT’s are used to obtain the Fourier representation from the 
data in physical space and thereby calculate derivatives. In order to avoid expensive evalua- 
tions of convolutions in Fourier space, the nonlinear terms are directly evaluated in physical 
space as products of derivatives. It has been shown by Canuto, Hussaini, Quarteroni and 
Zang (1988), that the ensuing aliasing error is negligible if the significant spatial scales of the 
computed variable are resolved on the grid. A third order, low storage Runge-Kutta scheme 
is used for advancing the solution in time. 

The computational mesh becomes progressively skewed because the mean flow velocity 
Ui with which the mesh moves has a variation in the x 2 direction. In order to control the 
loss of accuracy caused by excessive skewness of the mesh, we interpolate the solution on 
the grid at St = 0.5 shown in Fig. 2a, onto the grid sketched in Fig. 2b which is skewed 
in the opposite direction and corresponds to St = —0.5. Since the solution is periodic, the 
interpolation is straightforward in physical space. However, the remesh procedure introduces 
spectral aliasing errors as pointed out by Rogallo (1981). In the dealiasing procedure used 
here, the higher Fourier modes are truncated both before and after the remesh. 

Initial conditions have to be prescribed for u/, p : p and T. The initial velocity field is 
split into two independent components, that is, u/ = u \ f + uf \ each component having a 
zero average. The solenoidal velocity field u which satisfies V.u /J = 0 is chosen to be a 
random Gaussian field with the power spectrum 

E(k) — k A exp(— 2k 2 jk^) (5) 

where k m denotes the wave number corresponding to the peak of the power spectrum. The 
compressible velocity uf* which satisfies V X u c/ — 0 Is also chosen to be a random Gaussian 
field satisfying the same power spectrum, Eq. (5). The power spectra of the two velocity 
components are scaled so as to obtain a prescribed u TTnB = and a prescribed x = 
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u mia/ u rmj which is the compressible fraction of kinetic energy. The pressure p 1 * associated 
with the incompressible velocity is evaluated from the Poisson equation 


V p - -2pbu 2 l - pu { j u jti 


( 6 ) 


It remains to specify the thermodynamic variables. The mean density /? is chosen equal to 
unity, and p is chosen so as to obtain a prescribed Mach number 1 ^/ yj^ipj’p characterizing 
the turbulence. The fluctuating density p 1 and compressible pressure p c ' are chosen as 
random fields with the power spectrum Eq. (5) and prescribed pm* and p^. The pressure 
then becomes p = p + p 1 ' + p c \ the density is p — p + p', and the temperature T is obtained 
from the equation of state p = pRT. 


3 Results 

We have performed simulations for a variety of initial conditions and obtained turbulent fields 
with Taylor microscale Reynolds numbers Re\ up to 35 and turbulent Mach numbers M t up 
to 0.6. Note that Re\ = q\/v where q — and A = q/ , while M t = q/c where c 

is the mean speed of sound. The computational domain is a cube with side 2ir. The results 
discussed here were obtained with a uniform 96 3 mesh overlaying the computational domain. 
The initial values of some of the important non-dimensional parameters associated with the 
DNS cases ranged as follows: 3.6 < (SK/e ) o < 7.2, 16 < Re\ t0 < 24, 0.2 < M t< o < 0.4, and 
0 < x < 0.15. We note that K denotes the turbulent kinetic energy and e the turbulent 
dissipation rate. 

In previous work by Erlebacher et al. (1990) on isotropic turbulence, it was found that for 
a given initial Re \ , the initial choice of the non-dimensional quantities M t , Prm»/p, Tmu/T, 
Prma/p and % had a strong and lasting influence on the temporal evolution of the turbulence. 
The different choices of initial conditions that lead to different regimes were classified by 
Erlebacher et al. Even though the long time statistics in isotropic turbulence are, in general, 
strongly dependent on the initial conditions, there is a quantity in isotropic turbulence which 
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has a large degree of universality (in the absence of shocks). This quantity is the partition 
factor F defined by 

F = ~i 2 Mfx/Pc ( 7 ) 

which was shown by Sarkar et al. (1989) to equal unity (to lowest order) in a low M t 
asymptotic analysis, and to approach unity from a variety of initial conditions in direct 
simulations of isotropic turbulence with Mt up to 0.6. The physical interpretation of F —* 1 
is that there is a tendency towards equipartition between the kinetic energy and potential 
energy of the compressible component of the turbulence. 

In the case of initially isotropic turbulence subjected to homogeneous shear, for a given 
initial energy spectrum, the solution of the incompressible equations depends on the ini- 
tial conditions through the parameters ( SK/e) 0 and Re \ t 0 . The compressible case depends 
additionally on the initial levels of the thermodynamic fluctuations, dilatational velocity com- 
ponent and turbulent Mach number. Apart from the influence of initial conditions, another 
and perhaps more important feature is the influence of local compressibility and interactions 
between the velocity and thermodynamic fields on the statistics at a given time. Equilibrium 
scalings (if present) are also important to ascertain so as to improve the present capabilities 
of modeling compressible turbulence. In order to address these issues, we present selected 
results on second-order moments, pdf’s (probability density functions), and higher-order 
moments of both the velocity and thermodynamic fields. 

3.1 Second-order moments 

Evidence from physical and numerical experiments ( e.g., Tavoularis and Corrsin (1981) and 
Rogallo (1981) ) indicate that for incompressible homogeneous shear, both the turbulent 
kinetic energy and the turbulent dissipation rate increase exponentially. Though our sim- 
ulations are in accord with this picture of exponential growth, we note that Bernard and 
Speziale (1990) have recently proposed an alternative picture where, after a substantially 
long period of exponential growth, the turbulent energy eventually becomes bounded for 
large St through a saturation induced by vortex stretching. In order to explore the effect 


6 



of initial levels of compressibility on the evolution of the turbulent statistics, we have per- 
formed simulations based on two distinct types of initial conditions. The simulations of the 
first type, whose results are shown in Figs. 3a-6a, start with incompressible data, that is, 
Pnni ( o — Xo = 0? but have different initial Mach numbers M tt <> The second type of simula- 
tions, whose results are shown in Fig. 3b-6b, start with different levels of initial rms density 
\ 

fluctuations r Pt0 = (/3rms/p)o and compressible fraction of kinetic energy Xo> but have the 
same M t>0 = 0.3. We note that all the cases of Figs. 3-6 have Re\ t0 = 24 and (SK/e ) 0 = 7.2. 
The DNS results of Fig. 3a show that the level of the Favre-averaged kinetic energy K at 
a given time decreases with increasing M^o in the simulations with the first type of initial 
conditions, while Fig. 3b shows that the level of K also decreases with increasing %o and r Pf0 
for the simulations with the second type of initial conditions. Thus an increase in compress- 
ibility level, either due to increased Mach number or increased dilatational fraction of the 
velocity field decreases the growth of turbulent kinetic energy in the case of homogeneous 
shear flow. Fig. 4 shows that the development of the turbulent dissipation e has two phases. 
In the first phase (St < 4) the higher compressibility cases have higher values of e probably 
due to the buildup of the compressible dissipation e c . Later, in the second phase, the trend 
reverses and the compressible cases have lower e relative to the incompressible case. Thus, 
compressibility results in decreased growth of both e and K. 

In order to explain the phenomenon of reduced growth rate of kinetic energy, we consider 
the equation governing the kinetic energy of turbulence in homogeneous shear which is 

4(5 oK) = -S-pJfr'l -pe + ptf (8) 

at 

where ~pe is the turbulent dissipation rate and jfd' is the pressure-dilatation. The overbar 
over a variable denotes a conventional Reynolds average, while the overtilde denotes a Favre 
average. A single superscript ' represents fluctuations with respect to the Reynolds average, 
while a double superscript " signifies fluctuations with respect to the Favre average. It was 
shown in Sarkar et al. (1989) that the effect of compressibility on the dissipation rate can 
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be advantageously studied by using the decomposition 

e = e, + e c (9) 

where the solenoidal dissipation rate e, = uLu[ijj[ and the compressible dissipation rate e c = 
(4/3 )ud' 2 . Here uj[ denotes the fluctuating vorticity and d! denotes the fluctuating divergence 
of velocity. We note that the correlations involving the fluctuating viscosity have been 
neglected on the rhs of Eq. (9). Substituting Eq.(9) into Eq.( 8 ) gives 

^(P K ) = -Spu'lu'l - pe, - pe c + p 7 ^ (10) 

and the last two terms represent the explicit influence of the non-solenoidal nature of the 
fluctuating velocity field in the kinetic energy budget. The DNS time histories for e c and p'd' 
are presented in Figs. 5-6 respectively. The compressible dissipation e c is always positive (as 
expected) and when starting from a zero initial value as in Fig. 5a shows a monotonic increase 
with time. The pressure-dilatation jfd! is highly oscillatory as seen in Fig. 6 . Though p'd f 
can be both positive or negative it tends to be predominantly negative. The tendency of 
p*d f to be negative in homogeneous shear is clearly shown in Fig. 7 where the oscillations are 
smoothed out by plotting values of p'd f obtained by averaging over successive time intervals. 
Interestingly enough, on comparing Figs. 5a and 7a it seems that the time averaged values 
of —p 9 d* and e c are of the same order and exhibit similar trends of increase with St and M tf0 . 
These DNS results along with the kinetic energy equation ( 10 ) show that both e c and p f d f 
have a dissipative effect on the turbulent kinetic energy. 

The predominantly negative values of p'd* in homogeneous shear is in contrast to the 
predominantly positive values reported in Sarkar et al. (1989) for the case of decaying 
isotropic turbulence. To understand this contrasting behavior of p'd' we write the exact 
equation for the pressure variance which is applicable to both these flows 

= —2 7 p^ - (2 7 - l)pW- 2e p + <f> P (11) 

where <f> p is a term depending on mean viscosity and mean conductivity which is negligible 
compared to e p . Here e p = (7 — !)RKpTjTj and can be called the pressure dissipation 
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term by virtue of being a sink on the rhs of Eq. (11). Comparing Eq. (10) and Eq. (11) it 
is clear that p'd' acts to transfer energy between the kinetic energy of turbulence pK and 
the potential energy of the turbulence p' 2 /( 27p). In the case of homogeneous shear the rms 
values of the velocity and pressure increase with time. Now, the third term on the rhs of 
Eq. (11) is always negative and therefore a sink for the pressure variance. For small p' /p we 
can neglect the second term on the rhs of Eq. (11) with respect to the first term. Therefore 
for p' 2 to increase with time, a source term is necessary which implies that p'd' be negative 
in accord with the DNS results. We note that a similar analysis of the equations for p' 2 and 
T' 2 indicate that p'd 7 and T'd' also have to be predominantly negative in this flow. 

In decaying isotropic turbulence p' 2 has to decay with time; however, because e p is suffi- 
cient to ensure decay of p' 2 the sign of p'd' cannot be determined from Eq. (11). Alternatively, 
we consider the exact density variance equation which for homogeneous shear flow is 

— p* = -2pp^ - pycF (12) 

dt 

For small p'jp y the first term on the rhs of Eq. (11) dominates the second term. Therefore 
for p ' 2 to decrease with time it is necessary that p'd ' be predominantly positive. If we 
assume that the thermodynamic fluctuations are approximately isentropic in this case (DNS 
supports this assumption), it immediately follows that p'd' is also predominantly positive 
in decaying isotropic turbulence. Thus, the role of p'd' as a mechanism for energy transfer 
between the kinetic energy and potential energy dictates its differing signs in homogeneous 
shear and unforced isotropic turbulence. 

Fig. 8a shows the temporal evolution of Mt for different choices of the initial values for 
M t and Re A . For all the cases, the turbulent Mach number shows a monotonic increase 
with time for the range of St spanned by the simulations. Two results of the analysis of 
Sarkar et al. (1989) which were validated previously in isotropic turbulence are now checked 
in homogeneous shear flow. These results are, first, the compressible dissipation obeys the 
scaling e c /e 9 = aM f 2 , and second, the partition factor F ~ 1 (see Eq. (7) for the definition of 
F ) which implies a tendency toward equipartition between the kinetic and potential energies 
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of the fluctuating compressible mode. Fig. 9 a shows that for a variety of DNS cases the 
ratio e c /(e,M 2 ) approaches an equilibrium value of approximately 0.5 for large enough St. 
The same ratio is plotted as a function of local M t in Fig. 9 b. Thus Fig. 9 supports our 
previously derived result on the dependence of the compressible dissipation on turbulent 
Mach number. Zeman (1989) had proposed the scaling e c /e, = 1 — exp{[(M t — 0.1)/0.6] 2 } 
which is clearly in contradiction with the DNS results of Fig. 9. The partition factor F is 
seen in Fig. 10 a to approach and oscillate around an equilibrium value of approximately 0.95 
indicating that equipartition in the energies associated with the compressible mode holds in 
the case of homogeneous shear too. In Fig. 10 b we show the behavior of Ft = 7 2 M 2 x/p /3 
which is defined using the total pressure variance rather than the compressible pressure 
variance appearing in the definition of F. The parameter Ft also appears to reach values 
independent of initial conditions, however this trend is less pronounced in Ft relative to F. 

In Fig. 11 we show the time evolution of rms density and temperature fluctuations. It 
appears that both the rms normalized density and temperature fluctuations eventually scale 
as M 2 . The rms pressure also scales as M 2 . Erlebacher et al. (1990) decomposed the pressure 
p into an incompressible component p 1 and a compressible component p c . The behavior of 
p^ and is shown in Fig. 12. The instantaneous incompressible pressure p 1 ' satisfies 
the Poisson equation Eq. ( 6 ). The compressible pressure p c ' = p' — p 1 ' is the remaining 
part of the pressure which arises from both density variations and non-zero divergence of the 
velocity. Fig. 12 indicates that the rms of both components of the normalized pressure vary 
like Mf. A simple order of magnitude analysis of Eq. ( 6 ) along with recognition of the fact 
that SK/e = 0(1) shows that p^/p = 0(M 2 ). Thus it seems that in homogeneous shear 
the conventionally accepted Mach number variation holds for p^, p^ and Pmu* 

The evolution of p c ' in the case where the initial data is incompressible can be understood 
by considering the following equation for compressible pressure (see Erlebacher et al. (1990) 
for details) where heat conduction and viscous effects have been neglected 


V” + u!'r c :, + uf'p^ + uf' p r + 7(1 + P C 'K' = -( V + ufV'j) (13) 


10 



Since jF* and uf* are initially zero the growth of p c is due to the forcing by the incompressible 
pressure on the rhs of Eq. (13). This forcing causes p^^/p to become 0{^p TTaA l p) 0(M t ) on 
the convective time scale which explains the trend seen in the three lower curves of Fig. 12b. 
The top curve of Fig. 12b, which corresponds to a case with initially non-zero p Tms and 
also shows the same trend indicating that the forcing by the incompressible velocity field 
(through p 1 ) has a substantial influence in determining the level of 


3.2 Probability density functions and higher order moments 


The probability density function (pdf) gives a comprehensive description of a fluctuating 
variable because it formally enables the calculation of all the statistical moments. Second- 
order moments (discussed in the previous section) are probably the most important statistical 
indicators of the turbulence because they characterize intensities and cross-correlations of the 
fluctuating variables. Higher-order moments such as the skewness and flatness give additional 
information on the statistical distribution of the variable, for example, a positive skewness 
indicates a pdf which is skewed towards positive fluctuations from the mean, while a large 
flatness indicates a large propensity for the occurrence of extreme values of the variable. The 
skewness Sk(9) and flatness Fl(9) of a variable 9 are defined by 


Sk{9 ) = 
Fl(9) = 


9,3 

(py / 2 

P 

(Py 


A quantity with a Gaussian pdf has Sk = 0 and FI = 3. In this section, we present results on 
pdf, skewness and flatness of some variables for a DNS case with significant compressibility 
levels. The initial nondimensional parameters for this case are 


Rex = 24 , SK/e = 7.2 , Mt — 0.4 , X — 0 , Pma/ P — 0 j Prm*/P ® 


The pdf /„ of the streamwise fluctuating velocity in Fig. 13a appears to be symmetric 
around the origin. The skewness and flatness of u have insignificant deviations from Gaussian 
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values. The experiments of Tavoularis and Corrsin (1981) in homogeneous shear flow and 
Frenkiel and Klebanoff (1973) in the log-region of the flat plate turbulent boundary layer 
also indicate that u is approximately Gaussian. The pdf’s of the other components of 
the velocity are very similar to those of the streamwise velocity component. However the 
transverse velocity v develops a somewhat higher flatness ( Fl(y ) = 3.3 at St = 8 ) than 
the other velocity components. Fig. 13b shows the pdf f M of the Mach number M = qjc 
where q is the instantaneous flow speed and c is the instantaneous speed of sound. M has 
maximum values of 1.7 and 3.2 at St = 1 and St = 8 respectively, indicating that the pdf of 
M has exceptionally long tails. 

We now investigate quantities related to the velocity derivatives. Figs. 14a-b show the 
pdfs of the spanwise vorticity component cu z and the dilatation d. Both pdfs are skewed 
towards negative values. Higher-order moments of the vorticity components and dilatation 
are shown in Fig. 15. The skewness of the dilatation (see Fig. 15a) reaches values of about 
—0.5, signifying that local compressions are more likely than local expansions in the flow. 
Because of the invariance of the homogeneous shear problem to a reflection of the spanwise 
(z direction) axis, the cv x and u) y components have zero skewness. The dilatational field 
has larger flatness than the vorticity field, indicating that it is more intermittent than the 
vorticity field. 

In this DNS case, the maximum magnitude of dilatation in the compression regions 
attains values which are about an order of magnitude larger than the rms dilatation. Lee, 
Lele and Moin (1990) have also observed such high values of the maximum dilatation in 
isotropic turbulence, and have furthermore associated these high values with dynamically 
significant eddy shocklets. However, in order to assess the importance of these high dilatation 
regions it is necessary not only to identify the magnitude of the maximum dilatation but 
also to calculate the contribution of these regions to the rms dilatation. We find that, at 
St — 8, the compressive regions with a magnitude of dilatation larger than 3 times the rms 
dilatation contribute only about 2% to the variance of dilatation. This suggests that shocklet 
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dissipation is not the dominant contributory mechanism to the compressible dissipation for 

M t < 0.5. 

From Fig. 16a it appears that the imposition of shear does not alter the skewness and 
flatness of the Mach number M appreciably from the isotropic turbulence values at St = 0. 
The fact that M is bounded from below and not from above may be responsible for the 
departure of skewness and flatness from Gaussian values. The instantaneous values of M 
can become rather large ( the ratio = 5.8 at St = 8), but the flatness level of 4 

at St = 8 indicates that these large extrema have a relatively low probability of occurrence. 
The long tail in the pdf of M (see Fig. 13b and the previous discussion of the figure) does not 
seem to contribute significantly at the level of third and fourth moments. Fig. 16b shows that 
the velocity field can be highly supersonic ( M nuix = 3.3 at St = 8). Though the percentage 
of the volume occupied by supersonic fluctuations, denoted by <f> supeT i n Fig. 16b, increases 
with time, its maximum value is a modest 4.9% at St = 8. 

4 Characteristics of the instantaneous flow 

The previous sections have focused on the influence of compressibility on the statistical 
attributes of homogeneous shear turbulence. We now present some preliminary results on 
the spatial characteristics of the instantaneous turbulence from ongoing flow visualization. 
Figs. 17-20 are results from the flow field at St — 9 for a case which has the following initial 
nondimensional parameters 

Re x = 21 , SK/e = 5.8 , M t = 0.3 , x = 0.09 , p im ./p = 0.09 , p^/p = 0.13 

By St = 9 the initially isotropic turbulence develops into a realistic representation of ho- 
mogeneous shear turbulence. The turbulence at St = 9 has M t = 0.42, Prm*/p = 0.12 and 
X = 0.06, which implies a moderate level of compressibility, and has a Re\ = 32. One of 
the important characteristics of the velocity field is the extent and topology of supersonic 
regions. Fig. 17 shows a 3-D contour plot of regions with fluctuating Mach number M > 1. 
These regions of supersonic fluctuations are small discrete regions scattered through the flow 
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domain, and are elongated in the streamwise direction. The volume occupied by the regions 
with supersonic fluctuations Is small. Fig. 18 shows the dilatation V.u' in a plane parallel 
to the (xi,X 2 ) plane. The regions in black are expansion regions with positive dilatation, 
while the regions in grey and white are compression regions with negative dilatation. Both 
the compression regions and the expansion regions are very elongated in the streamwise x 1 
direction. The white portions in the dilatation plot of Fig. 18 correspond to compression 
regions where the magnitude of dilatation is larger than the maximum positive dilatation. 
The magnitude of dilatation in these white regions are up to 50% larger than the maximum 
positive dilatation indicating that compressions are stronger than the expansions. As Fig. 18 
shows, the white regions of strong compression are small and discrete elements which occupy 
a small portion of the plane. Figs. 19a and 19b are contour plots of the magnitude of u and 
uf respectively in a plane parallel to the (xi,x 2 ) plane. Regions with large gradients seem 
to be more prevalent in the uf * field relative to the u field. Figs. 20a and 20b show the 
incompressible pressure p 1 ' and the compressible pressure p c * in the same plane considered 
in Fig. 19. It is clear that the p 1 ' field has significant differences with respect to the p c * field. 
The compressible pressure has higher gradients and more pronounced small scale features 
than the incompressible pressure. Elongated and inclined structures are also visible in the 
p c f field. 

5 Conclusions 

We have performed direct numerical simulations of homogeneous shear flow using a spectral 
collocation technique and analyzed the data base in order to quantify compressibility effects 
on the turbulence. An increase in the compressibility level (i.e. turbulent Mach number, 
density fluctuations or rms dilatational velocity fluctuations) leads to a decrease in the 
growth of turbulent kinetic energy. In the compressible case, the kinetic energy equation 
has two additional terms, the compressible dissipation and the pressure-dilatation. Both 
these terms have a dissipative contribution in shear flow, leading to the reduced growth of 
the turbulent kinetic energy. The O(M^) variation of the compressible dissipation which 
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was obtained previously by an asymptotic theory and then verified in the DNS of isotropic 
turbulence, has now been substantiated by the DNS of homogeneous shear turbulence. The 
pressure-dilatation is a more difficult term to characterize, because it can be either positive 
or negative and also because gradient transport arguments are clearly inapplicable. Our DNS 
show that the pressure-dilatation is predominantly negative in the homogeneous shear case 
in contrast to its predominantly positive character in decaying isotropic turbulence. Analysis 
shows that, in the transport equation for the pressure variance, the pressure-dilatation has 
to be a source in the case of homogeneous shear (where the turbulence grows) and a sink 
in the case of decaying isotropic turbulence, thus accounting for the opposite signs of the 
pressure-dilatation in the two flows. In the core region of inhomogeneous aerodynamic flows, 
perhaps the pressure-dilatation is predominantly negative so as to transfer energy from the 
velocity field to the pressure field. 

It is well known that, for low-speed flows, the rms pressure normalized by the mean 
pressure varies as 0(A/ t 2 ). Our DNS results support this scaling of the rms of pressure, 
density and temperature fluctuations for compressible homogeneous shear flow. These DNS 
results and our analysis suggest that such a scaling may be valid for those inhomogeneous 
flows where the thermodynamic fluctuations arise due to high-speed effects and not due 
to other factors such as external heating or cooling, and mixing of fluids with different 
densities or temperatures. We have studied the probability density functions and higher- 
order moments of the dilatation and the vorticity. The dilatation field has appreciable 
negative skewness indicating a preponderance of compressions over expansions, and has a 
larger flatness than the vorticity indicating that it is more intermittent than the vorticity. 
The intermittent regions of negative dilatation were found to have only a small contribution 
to the compressible dissipation for M t < 0.5. 
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Figure 3 Evolution of kinetic energy for various DNS cases. 
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Figure 6 Evolution of pressure-dilatation for various DNS cases 
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Figure 14a Pdf of spanwise vorticity component. 
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Figure 14b Pdf of dilatation. 
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Figure 15a Skewness of vorticity components and dilatation. 
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Figure 15b Flatness of vorticity components and dilatation. 
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Figure 16 a Skewness and flatness of Mach number. 



St 


Figure 16 b Maximum Mach number and volume fraction of supersonic fluctuations. 
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Figure 1 7 Perspective of supersonic regions in the flow. 
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Figure 19a Contours of the magnitude of u\' in a (xijXa) cut of the flow domain. 
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Figure 19b Contours of the magnitude of uf' in the plane shown in Fig. 19a. 
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Figure 20a Contours of p v in the plane shown in Fig. 19a. 
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